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Abstract 

We show that in brane- world scenarios with warped extra dimensions, 
the Casimir force due to bulk matter fields may be sufficient to stabilize the 
radion field (p. In particular, we calculate one loop effective potential for 
(j) induced by bulk gravitons and other possible massless bulk fields in the 
Randall-Sundrum background. This potential has a local extremum, which 
can be a maximum or a minimum depending on the detailed bulk matter 
content. If the parameters of the background are chosen so that the hierarchy 
problem is solved geometrically, then the radion mass induced by Casimir 
corrections is hierarchycally smaller than the TeV . Hence, in this important 
case, we must invoke an alternative mechanism (classical or nonperturbative) 
which gives the radion a sizable mass, to make it compatible with observations. 
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I. INTRODUCTION 



It has been suggested that theories with extra dimensions may provide a solution to 
the hierarchy problem The idea is to introduce a dimensional internal space of 

large physical volume V, so that the the effective lower dimensional Planck mass rripi ~ 
yi/2j^^(rf+2)/2 -g j2iuch larger than M ~ TeV- the true fundamental scale of the theory. In 
the original scenarios, only gravity was allowed to propagate in the higher dimensional bulk, 
whereas all other matter fields were confined to live on a lower dimensional brane. Randall 
and Sundrum (RS) introduced a particularly attractive model where the gravitational 
field created by the branes is taken into account. Their background solution consists of two 
parallel flat branes, one with positive tension and another with negative tension embedded 
in a a five-dimensional Anti-de Sitter (AdS) bulk. In this model, the hierarchy problem is 
solved if the distance between branes is about 37 times the AdS radius and we live on the 
negative tension brane. More recently, scenarios where additional fields propagate in the 
bulk have been considered 0-0. 

In principle, the distance between branes is a massless degree of freedom, the radion 
field 0. However, in order to make the theory compatible with observations this radion must 
be stabilized [|^-[T0[|. Clearly, all fields which propagate in the bulk will give Casimir-type 
contributions to the vacuum energy, and it seems natural to investigate whether these could 
provide the stabilizing force which is needed. In this paper, we shall calculate the radion 
effective potential V^jj^cj)) due to KK gravitons and other massless bulk fields. As we shall see, 
this effective potential has a rather non-trivial behaviour, which generically develops a local 
extremum. Depending on the detailed matter content, the extremum could be a maximum 
or a minimum, where the radion could sit. For the purposes of illustration, here we shall 
concentrate on the background geometry discussed by Randall and Sundrum, although our 
methods are also applicable to other geometries, such as the one introduced by Ovrut et al. 



in the context of eleven dimensional supergravity with one large extra dimension fTT 



The paper is organized as follows. In Section II, we introduce the model and conventions. 
In Section III we calculate the effective potential due to conformally invariant bulk fields, 
and discuss the stabilization mechanism. Conformally invariant fields are convenient because 
of their simplicity, which allows a quick derivation of their contribution. Moreover, the 
backreaction of the Casimir energy on the geometry can be taken into consideration in this 
case. Section IV is devoted to the graviton contribution. For gravitons, conformal invariance 
is lost, and the evaluation of the effective potential is not straightforward. A method will be 
developed which relates the curved space effective potential to a suitable effective potential 
in flat space. The later takes the form of a sum over the contributions of an infinite tower 
of massive KK fields, whose mass spectrum m„(0) is a function of the brane separation (the 
discrete index n labels the infinite KK tower). Our conclusions are summarized in Section 
V. 

Related calculations of the Casimir interaction amongst branes have been presented in an 
interesting paper by Fabinger and Horava |]T2|. In the concluding section we shall comment 
on the differences between their results and ours. 
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II. THE RANDALL-SUNDRUM MODEL AND THE RADION FIELD 



To be definite, we shall focus attention on the brane-world model introduced by Randall 
and Sundrum [Q. In this model the metric in the bulk is anti-de Sitter space (AdS), whose 
(Euclidean) line element is given by 

ds^ = o?{z)riabdx''dx^ = o?{z) dz"^ + dx^ = dy^ + {z)dx^ . (2.1) 

Here a{z) = i/z, where £ is the AdS radius. The branes are placed at arbitrary locations 
which we shall denote by 2+ and z^, where the positive and negative signs refer to the 
positive and negative tension branes respectively {z^ < z_). The "canonically normalized" 
radion modulus - whose kinetic term contribution to the dimensionally reduced action on 
the positive tension brane is given by 

\j dS^g%'d,(t)d,<j), (2.2) 

-is related to the proper distance d = Ay between both branes in the following way 

= (3M3£/47r)i/'e-'^/^ 

Here, M ~ TeV is the fundamental five- dimensional Planck mass. It is usually assumed 
that i ~ . Let us introduce the dimensionless radion 



\3MHJ 



Z- 



which will also be refered to as the hierarchy. The effective four-dimensional Planck mass 
rupi from the point of view of the negative tension brane is given by rUpi = M^i{X~^ — 1). 
With d ~ 37i, A is the small number responsible for the discrepancy between nipi and M. 

At the classical level, the radion is massless. However, as we shall see, bulk fields give rise 
to a Casimir energy which depends on the interbrane separation. This induces an effective 
potential X;jf{<P) which by convention we take to be the energy density per unit physical 
volume on the positive tension brane, as a function of 0. This potential must be added to 



the kinetic term ( p.2|) in order to obtain the effective action for the radion: 

■1 



Sgff [4>] = I d xttj 



(2.3) 



In the following Sections, we calculate the contributions to V^jj from conformally invariant 
bulk fields and from bulk gravitons. 



III. RADION STABILIZATION BY MASSLESS BULK FIELDS 

Let us start by considering the contribution to \^jj{(f)) from conformally invariant mass- 
less bulk fields. Technically, this is much simpler than finding the contribution from bulk 
gravitons which will be discussed in the next section. Moreover the problem of backreaction 
of the Casimir energy onto the background can be taken into consideration in this case. 
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Consider, for instance, a conformally coupled scalar x- This obeys the equation of motion 

-□,X + ^^^i?X = 0. (3.1) 



where is the d'Alembertian operator in the metric and R is the Ricci scalar. Here, 
we consider the case of arbitrary odd spacetime dimension Z), with branes of co-dimension 
one. By changing the variable X X = c'-'^^^^-'^^X, the equation of motion for x becomes 

□ X = 0. (3.2) 

Here □ is the flat space d'Alembertian. It is customary to impose Z2 symmetry on the bulk 
fields, which results in Neumann boundary conditions 

d.X = 0, 

at z+ and z_. The eigenvalues of the d'Alembertian subject to these conditions are given by 

>l,k={^y + k\ (3.3) 

where n is a positive integer, L = — z+ is the coordinate distance between both branes 
and k is the coordinate momentum parallel to the branes. 

Similarly, we could consider the case of massless fermions in the RS background. The 
Dirac equation is conformally invariant ||13| , and the conformally rescaled components of the 
fermion obey the flat space equation ( ^.2[ ) with Neumann boundary conditions. Thus, the 
spectrum ( p.3|) is also valid for massless fermions. 

Let us now consider the Casimir energy density in the conformally related flat space 
problem. We shall first look at the effective potential per unit area on the brane, A. For 
bosons, this is given 

= ^Tr ln(-n^^). (3.4) 

Here /i is an arbitrary renormalization scale. Using zeta function regularization [See e.g. 
I^, or the discussion in Section IV for details], it is straightforward to show that 



Here t] = (D — l)/2, and Qr is the standard Riemann's zeta function. The contribution of 
a massless fermion is given by the same expression but with opposite sign: 

Vi{L) = -V^{L). (3.6) 

The expectation value of the energy momentum tensor is traceless in flat space for confor- 
mally invariant fields Moreover, because of the symmetries of our background, it must have 
the form |]TE| 

{T\)^lat = {D- l)po(^), {Ty)siat = -Poiz) 5V 
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By the conservation of energy-momentum, po must be a constant, given by 



^° 2L ^2L«' 

where the minus and plus signs refer to bosons and fermions respectively and we have 
introduced 

(47r)''r/! ' 

Now, let us consider the curved space case. Since the bulk is of odd dimension, there is 
no anomaly [|T^] and the energy momentum tensor is traceless in the curved case too. This 
tensor is related to the fiat space one by (see e.g. |[T3[] ) 



Hence, the energy density is given by 

p = a~^po. (3.7) 

The effective potential per unit physical volume on the positive tension brane is thus given 
by 

V,ff{X) = 2 aX-^J a''{z)pdz = T^^-^ _ . (3.8) 

Note that the background solution a{z) = i/z has only been used in the very last step. 

The previous expression for the effective potential takes into account the casimir energy of 
the bulk, but it is not complete because in general the effective potential receives additional 
contributions from both branes. As we shall discuss in more detail in Section IV, we can 
always add to V^jj terms which correspond to finite renormalization of the tension on both 
branes. These are proportional to A° and A^~^. The coefficients in front of these two 
powers of A cannot be determined from our calculation and can only be fixed by imposing 
suitable renormalization conditions which relate them to observables. Adding those terms 
and particularizing to the case of = 5, we have 



%(A) = Tt 



[i-xy 



(3.9) 



where A ^ 2.46 ■ 10"'^. The values a and (3 can be obtained from the observed value of 
the "hierarchy", Xobs, and the observed value of the effective four-dimensional cosmo logical 
constant, which we take to be zero. Thus, we take as our renormalization conditions 

VeffiXobs) = -^i^obs) = 0. (3.10) 

If there are other bulk fields, such as the graviton, which give additional classical or quantum 
mechanical contributions to the radion potential, then those should be included in V^jj. From 
the renormalization conditions ( |3.10| ) the unknown coefficients a and (3 can be found, and 
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FIG. 1. Contribution to the radion effective potential from a massless bulk fermion. This is 
plotted as a function of the dimensionless radion A = e'^'^^^ , where d is the physical interbrane 
distance. The renormalization conditions ( ^.10 ) have been imposed in order to determine the 
coefficients a and (3 which appear in (f 



then the mass of the radion is calculable. In Fig. 1 we plot (|3.9|) for a fermionic field and a 
chosen value of \obs- 
From ( p.lO| ), we have 



(3 = -A{l-X 



obs I 



a 



obs' 



(3.11) 



These values correspond to changes 6a±. on the positive and negative brane tensions which 
are related by the equation 



(3.12) 



As we shall see below, Eq. ( |3.12| ) is just what is needed in order to have a static solution 
according to the five dimensional equations of motion, once the casimir energy is included. 

We can now calculate the mass of the radion field m^^^ from the point of view of the 
negative tension brane. For Xobs ^ 1 we have: 



m 



2 (-) 



A 



obs "^(j, 



2 (+) 



A 



obs 



obs 



5vr^a(-4) 
6M3/5 



(3.13) 



The contribution to the radion mass squared is negative for bosons and positive for fermions. 
Thus, depending on the matter content of the bulk, it is clear that the radion may be 
stabilized due to this effect. 
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Note, however, that if the "observed" interbrane separation is large, then the induced 
mass is small. So if we try to solve the hierarchy problem geometrically with a large internal 
volume, then Xobs is of order TeV/rUpi and the mass ( p.l3| ) is much smaller than the TeV 



scale. Such a light radion would seem to be in conflict with observations. In this case we 
must accept the existence of another stabilization mechanism (perhaps classical or nonper- 
turbative) contributing a large mass to the radion. Of course, another possibility is to have 
Xobs of order one, with M and i of order rripi, in which case the radion mass would be 

very large, but then we must look for a different solution to the hierarchy problem. 

Due to conformal invariance, it is straightforward to take into account the backreaction of 
the Casimir energy on the geometry. First of all, we note that the metric is analogous to 
a Friedmann-Robertson- Walker metric, where the nontrivial direction is space-like instead of 
timelike. The dependence of a on the transverse direction can be found from the Friedmann 
equation 

fa'Y levrGg A 



Here a prime indicates derivative with respect to the proper coordinate y [see eq. (|2.1| )], 
and A < is the background cosmological constant. Combined with ( p.7| ), which relates the 
energy density p to the scale factor a, Eq. (|3.14|) becomes a first order ordinary differential 
equation for a. We should also take into account the matching conditions at the boundaries 

f a'\ =F87rG'5 
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-a±. (3.15) 



A static solution of Eqs. (|3.14|) and (|3.15| ) can be found by a suitable adjustment of the 



brane tensions. Indeed, since the branes are fiat, the value of the scale factor on the positive 
tension brane is conventional and we may take a+ = 1. Now, the tension a+ can be chosen 
quite arbitrarily. Once this is done, Eq. ( |3.15| ) determines the derivative a^, and Eq. ( ^.141 ) 



determines the value of po- In turn, pq determines the co-moving interbrane distance L, 
and hence the location of the second brane. Finally, integrating ( p.l4| ) up to the second 
brane, the tension a- must be adjusted so that the matching condition ( |3.15[ ) is satisfied. 



Thus, as with other stabilization scenarios, a single fine-tuning is needed in order to obtain 
a vanishing four-dimensional cosmological constant. 

This is in fact the dynamics underlying our choice of renormalization conditions (|3.10|) 



which we used in order to determine a and [3. Indeed, let us write a+ = ao + 5cr+ and 
(T_ = —(To + 6a-, where ctq = (S/iniG^) is the absolute value of the tension of the branes in 
the zeroth order background solution. Elliminating a' /a from ( p.l5|) and ( p.l4|) , we easily 



recover the relation ( p.l2|) , which had previously been obtained by extremizing the effective 
potential and imposing zero effective four-dimensional cosmological constant (here, 6a± is 
treated as a small parameter, so that extremization of the effective action coincides with 
extremization of the effective potential on the background solution.) In that picture, the 
necessity of a single fine tuning is seen as follows. The tension on one of the walls can be 
chosen quite arbitrarily. For instance, we may freely pick a value for (3, which renormalizes 
the tension of the brane located at . Once this is given, the value of the interbrane distance 
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Xobs is fixed by the first of Eqs. Then, the value of a, which renormahzes the tension 

of the brane at z^, must be fine-tuned to satisfy the second of Eqs. (|3.1CI|) . 

Eqs. ( p.l4| ) and ( |3.15| ) can of course be solved nonperturbatively. We may consider, for 
instance, the situation where there is no background cosmological constant (A = 0). In this 
case we easily obtain 



6nGA 

{z_ - z+y 



(3.16) 



where the brane tensions are given by 2nGa± = ±(C — z±)~^ and C is a constant. This is a 
self-consistent solution where the warp in the extra dimension is entirely due to the Casimir 
energy. 



IV. GRAVITON CONTRIBUTION TO THE RADION EFFECTIVE POTENTIAL 

Since the graviton is not conformally invariant, the techniques needed in order to find 
its contribution to the effective potential will have to be somewhat more sophisticated. 
However, because of the relevance of the graviton, it seems worthwile to pursue this task 
and in the following subsections we develop a convenient method to address this calculation. 
This consists of two steps. In Subsection IV.A we relate the curved space effective potential 
to a suitable fiat space determinant. The calculation of this determinant is then presented 
in Section IV.B. 

Our background has maximally symmetric foliations orthogonal to the y direction, and 
just like in the case of cosmological Friedmann-Robertson- Walker models, each graviton 
polarization contributes as a massless minimally coupled scalar field. Although this seems 
to be a well known fact which is often quoted in the literature, it is not always clear from 
the context what is the extent of this equivalence. The correspondence is straightforward at 
the classical level. At the quantum level, it can also be shown to hold, although this is not 
so straightforward to prove because careful gauge-fixing of the gravitational sector must be 
done. A detailed justification is lengthy, and will be presented elsewhere |I5| . 

In short, the correspondence is as follows. Perturbations of the gravitational field in the 
Randall-Sundrum gauge can be expressed as 

/^a6 = ^?a6-^7i?=E«M?%), (4-1) 

where CT55 = a^i, = and o"^,^ is a constant transverse and traceless polarization tensor 
(a Fourier decomposition of $ is assumed along the branes, in order to define transverse 
polarization). The summation is taken over all polarizations. The quadratic reduced action 
for one particular polarization becomes 

W = J d^x^/g^ {-^g^) + (boundary term), (4.2) 
where is the usual covariant scalar Laplacian associated with the five dimensional metric 

Qab = a^Vab, 
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and we omitted the index (i). 

Thus, we shall be interested in the determinant of the operator 



P 



-□ 



5' 



(4.3) 



with appropriate boundary conditions at the branes. These are determined from Israel's 
junction conditions plus the requirement of Z2 symmetry |TB[. Equivalently, they can be 



found from the equation of motion for $ in the RS background. It is easy to show that in 
terms of $, they reduce to the standard Neumann boundary conditions i.e., 



9.$ = 0, 



(4.4) 



at z = z±. 



A. From flat to curved space 

Our next task is to evaluate the determinant of the covariant d'Alembertian operator P 
defined in (|4.3|). Formally, this is the product of its eigenvalues. Nevertheless, in the present 



case, a straightforward calculation involving the eigenvalues of P seems rather complicated 
and not particularly illuminating. For instance, the eigenvalues of P are completely unrelated 
to the spectrum of Kaluza-Klein gravitons, and so the contribution from individual physical 
modes seems hard to identify. 

A more promising strategy is to express (det P) in terms of the determinant of an operator 
whose eigenvalues are directly related to the KK masses. For this purpose, we introduce 
a one parameter family of metrics which interpolates between flat space and our physical 
space. 

For definiteness we shall take 

where the parameter e runs from to 1. For e = the fictitious metric cjab is fiat space, 
whereas for e = 1 it coincides with the physical metric . The arguments presented in this 
section are independent of the precise form of the "path" (see Fig. 2). The explicit form 
( [4.5[ ) will only be used at the very end of the calculation, and it is chosen for convenience 
so that all metrics in the path are AdS spaces with curvature radius i/e. 
Next, we introduce the operator P^ defined by 

nr^p^ni^ = nfp = (^z'^-'d^z'-^d, + di^ . (4.6) 

This operator can be written in covariant form as 

P, = -(a, + E^, 
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e=0 FLAT 

FIG. 2. Paths from flat to AdS space in the space of conformal factors. The mtegral in ( [4.15| ) is 
independent on the path. 

where is the d'Alembertian in the space with metric Qab and 

2^ + {D-4)g'^''{dJnn,){d,\nn,) . 
For later reference, we note that acts on the rescaled field 

<i. = a^$, (4.7) 
which obeys the boundary condition B^^ = {n°-da + Sf)^\dM = 0. Here 

S, = ^^n''djnn,, (4.8) 

and n is the inward directed normal vector at the boundary. 

We are interested in the effective potential per unit co- moving area Ac- This is given by 
= A^^ ln(det Py/"^. As we shall see, the dependence of (det P^) on e along the conformal 
path in Fig. 2 is easily found in terms of local quantities. Hence, instead of dealing directly 
with (detP), we shall calculate the fiat space determinant (detPo); and then integrate the 
dependence on e along the path. The effective potential per unit comoving area is given by 

Veff = 

where 



V 



4fi2 
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Here, we have introduced 



-Tr In 



1 



(4.9) 



C(z/) = Tr 



2/i 



2;^ 



(4.10) 



The symbol Tr refers to the usual L2 trace. For any operator (9, the trace can be represented 
as 

Tr[0] = Y^jd'>x g^'^ $,(0$,) = T.jd''x f'^ <l>.(0<l.), 

i i 

where $j (or $j) form an orthonormal basis with respect to the measure associated with the 
metric g (or g respectively). It is straightforward to show that 



Tr 



Tr 



(4.11) 



Following the Ref. |17-20|, we introduce the asymptotic expansion of the trace for small 



Tr 



/(x)e 



Er-V/2(/,Pe), 



(4.12) 



n=0 



where / is an arbitrary smooth function. This is a generalization of the more widely used De 
Witt-Schwinger expansion for the solution of the heat kernel equation |^ (which corresponds 
to the case / = 1). The dependence of K on e is related to the coefficient a£)/2 in the following 



way 



9,lim9X.(^) = limS,-^ T c^^'^Tr (-/.S^e-^^ 



lim 



- r(z/) 

lim d,y - 



-2U-D-1 



-ePe 



-0 "{D- 2u){D - 1 - 2z/) ■ ■ ■ (1 - 2z/; 
2a^/2(/„P.), 



-'^ n=0 



(4.13) 



where we have defined 



(4.14) 



In the first equality, we have used ( [4.11| ). In the second equality, we performed integration 
by parts assuming 2v > D. The contribution from ^ = cxd vanishes due to the exponential 
suppression in Tr f{x)e~^ . In the third equality, we performed a further integration by 
parts. 

The dependence in e can now be integrated along the path which joins fiat space with 
AdS, and the effective potential can be expressed as 
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V^^ = V^ = Vo-^ J^'deaD/2{fe,Pe). 



(4.15) 



The calculation of Vq is done in the next section. The coefficient ci_d/2, for D < 5, has 
been studied in |T^J20J2^-p5| for a large class of "covariant" operators which includes our 



Pg. In general, this can be expressed in terms of integrals over the bulk and boundary of 
local invariants constructed from the metric g and the functions 5*^ and E^^ which enter the 



boundary condition, given in Eq. (^ 



For the case of odd dimension D, the coefficients ac/2 contain only contributions from 
the boundary. For instance, for D = 3 we have p2 



a3/2(/, Pe) = f {/(96E, + IQR - 8P^^„ + 13K^ + 2K^,K^'' 

153671 JdM 

-9QS,K + 192^2) + /;„(-6K + 9QS) + 24/.„„}. 

Here K^^, is the extrinsic curvature and i?^^^ = +Thdc ~ ■■■ the Riemann tensor. As usual, 
greek indices run from to 3 and the index n indicates contraction with the inward unit 
normal vector h. If we integrate this expression along the path given by ( [1.15| ), then all 
intermediate geometries are AdS of radius i/e and the invariants are easily found. The 
result is given by 



j-/'c^ea3/,(/.,P.: 



1 



l(lnz+ + AMn._)-A(i + A2) 



(4.16) 



where we have reintroduced the dimensionless radion A 
for 05/2, relevant to the five dimensional case, has been given by Kirsten |26 



The general expression 
It is rather 



lengthy since it involves combinations of geometric invariants that contain 4 derivatives of 
the metric, and we shall not reproduce it here. Integrating this coefficient over e, we obtain 



-l/^ea5/2(/.,P. 



2^4 



27 53 
128^ + ^ 5120^ 



(4.17) 



It should be noted that under the rescaling z± — ^ -t^z^ (with ^ = const.), which corresponds 
to a translation of both branes in the y direction, the expression for the effective potential 
V^jj must scale like 



V%: 

eff 



^1 



eff 



(4.18) 



This shift is a mere coordinate transformation which changes the physical volume corre- 
sponding to a given co-moving volume. We shall show in the next section that the quantity 
Vq does not scale in this way, and that the integrals ( [4. 161) and (|4.17|) have precisely the form 
which restores the scaling ( |4.18| ). For that reason, the explicit evaluation of these integrals 
is not strictly necessary, although of course it provides a valuable consistency check. As we 
shall see, the integrals can be guessed from the scaling property ( 4.18| ) up to terms which 
can be reabsorbed by finite renormalization. 
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B. Computation of %ff{(p) 



From ( [4.6|) , the eigenvalue problem for e = takes the form: 



(4.19) 



with boundary condition dz^i = at z = z±. Here, an expansion in Fourier modes of 
wavenumber k along the flat branes has been assumed. In terms of the rescaled field $ defined 
this equation looks like a Schrodinger equation with a "volcano" -type potential |]16 



m 



Defining mf = — fc^ + £^Aj, we have 



.D-2; 



= 0, 



(4.20) 



which has a general solution $j = z"^ (Ai J^(mjz) + A2Yrj{miz)) . The boundary conditions 
give an equation which determines rrii as 



F{mi) = J^_i(miA)F^_i(mi) - y;,_i(miA) J^_i(mj) = 0, 



(4.21) 



where J and Y are Bessel functions, rj = {D — l)/2, rhi = miZ_ and we have reintroduced 

A = z+/z_. 

Equation ( [4.21j ) gives the physical spectrum of masses m_ = rhi/i for the KK gravitons 
from the point of view of the negative tension brane. This spectrum is a function of the 
interbrane distance, which enters Eq. ( ^.21| ) through the hierarchy A. The zeta function for 
the operator Pq, given by 



Co(^)=A/ 



d 



D-lj 



D 



(4.22) 



resembles the zeta function for a tower of fields of mass m-j in flat space. Nevertheless, 
this is to some extent just a formal analogy: instead of physical momentum and mass, the 
coordinate momenta k and "coordinate" mass eigenvalue mi appear in this expression. Q 
Performing the momentum integrals in ( 4.22|) we have 



Co(^^) = A 



{ifiY'^z 



T{v-ri) 



(47r)'?r(z/) 



C(2z/-(D-1)), 



(4.23) 



^In fact, it does not seem possible to write down an expression analogous to (4.22) for our V^g, 
containing physical masses and momenta instead of coordinate ones. This would be possible if our 
space-time were ultrastatic (see e.g. [^]), which just means that the lapse function is a constant. 
However, in our case, time flows at different rates in different places, so for instance the physical 
mass depends on which brane we are considering. This seemingly trivial complication prevents us 
from doing the naive "projection" of the problem onto four dimensional slices, and renormalizing 
infinities as if we were considering an infinite tower of fields living in a fiat space. In fact, such 
naive procedure would give a qualitatively different result for X;jfi^) depending on whether we 
compute it on the positive or on the negative tension brane, which is clearly nonsensical. This is 
why the techniques of Section IV. A have been introduced. 
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FIG. 3. Contour of integration in Eq. (|4.26| ). 



where we defined 



c(^) = E 



(4.24) 



Substituting ( ^^231) in ([4.9|) we have 



Vo{X) 



D-l 



I 



^ln(£/i^_) + 2 E [ C(l -D) + C'(l - D) 



n=l 



(4.25) 



Now the problem reduces to the calculation of the generalized zeta function ( [4.24| ). This is 
a discrete sum over the zeros of the function F defined in (|4.21| ). This type of sums, where 
F is a combination of Bessel functions, have been previously studied in [p^ -|3T|l. Following 
29| , we express ({s) as 



1 r F' H r 

ns) = — / dtr'— = — / dtt- 

' 2m Jc F 2m Jc 



l-s 



InF. 



(4.26) 



This follows from the fact that F has only simple zeros and these are on the positive real 
axis. The contour of integration C is given in Fig. 3. The expression (|4.26| ) is convergent 
for sufficiently large s. 

After some manipulations, which are deferred to Appendix A, we find 



Vo{X) 



D-l 



Ik + X^^-'li + /3d-i(1 + A^-i) | ln(£/i) + ^ J] n'^ 

^ n=l 
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FIG. 4. The function V(A). 



+I3d-i (in z+ + A^-i In z_) + A^^^ V(A) 
Here Tj^Xk and /3z)_i are constants which are evaluated in the Appendix, and 



(4.27) 



where / and K are the modified Bessel functions. The function V(A) is plotted in Fig. 4. 

The coefficients which accompany the powers A° and A"^"^ will change if we change the 
renormalization scale /z inside the logarithm. This corresponds to the well known scaling of 
the effective potential under changes of /i, which adds a geometric term proportional to the 
coefficient aD/2{f = 1,-P) defined in (|4.16|) . It is easy to show that this coefficient is Weyl 
invariant and therefore independent of e. The invariant is made out of local terms (such as 
the square of the extrinsic curvature) which should also be included in the bare action. The 
renormalized couplings in front of these local terms are in fact supposed to be negligibly 
small (they have not been included in the equations of motion satisfied by the background), 
so the terms proportional to In fi might as well be dropped from Vq. On the other hand, we 
can always add to Vq terms which correspond to a finite renormalization of the cosmological 
constant in the bulk. These will give contributions proportional to (1 — A"^"^). Also, one can 
add terms which correspond to a finite renormalizaton of the tension on the branes. These 
will give independent contributions proportional to A° from the positive tension brane and 
to A^^^ from the negative tension brane. Hence, the coefficients in front of these two powers 
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of A cannot be determined from the calculation and can only be fixed by imposing suitable 
renormalization conditions which relate them to observables. 

By using ( |4.15| ) with the results ( [4.16| ) or ( |4.17| ) obtained in the preceding section, we can 
now evaluate V^'jj. As mentioned at the end of the last subsection, the explicit calculation 
of the generalized Seeley-De Witt coefficients is not strictly necessary. It is known that aD/2 
for odd D has only contribution from the terms depending on the background quantities 
evaluated on the boundaries. Therefore, the shift in the effective potential due to the 
variation of e must be given by the sum of a function of plus a function of z_: 



1 /•! 

— deaD/2{fe,P,) = ^+{z+) + T-{zJ). 



(4.29) 



Imposing that V^^ has the correct scaling behaviour given in Eq. (|4.18|) , the functional form 

J-± is determined up to the two coefficients in front of A° and A"^"^ mentioned above, which 
have the same A dependence as cosmological constants on the respective branes. Taking this 
into account, we have 



Vie 



(47r)^(r/ - iy.z_ 



D-l 



a + /3A^-^ + A^-^V(A) 



(4.30) 



The "renormalized" values a and f3 can be obtained as before from the renormalization 
conditions (|3.10|). An example of an effective potential after these conditions have been 



imposed is shown in Fig. 5. Here 

%(A)=ai-^(^+)l^^(A), 

is the effective potential per unit physical volume on the positive tension brane, which 
appears in Eq. ( p.3|) . (Note that, in fact, the conditions ( p.lO|) do not depend on wether we 
are using the effective potential per unit co- moving or per unit physical volume.) 
Expanding V(A) for small A and in D = 5, we have 



V(A) = JA2 + 0(AMnA) 



where 



I=--l dpp'^^'^P^ 



1 /-^o 

2 Jo 



hip) 



-4.2. 



(4.31) 



(4.32) 



Hence, 



%(A) 



(Any 



(4.33) 



The second equality in (|3.1CI| ) determines f3 ~ —{3/2)IXlf^g, and then the first of the renor- 
malization conditions gives the aditive constant a. With this, the physical mass of the 
radion from the point of view of the negative tension brane is given by 



m 



2 (-) 



^obs 



TC 



I 



(4.34) 
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FIG. 5. Graviton contribution to the radion effective potential. 



For the case of the graviton considered here, this mass squared is negative (this is shown 
here for small Xots, but it is in fact true for all Xobs), and therefore the graviton cannot by 
itself stabilize the radion. 

As we have seen in Section II, in theories where additional fields can propagate in the 
bulk, some of them may give positive contributions to the mass squared. In fact, for a small 
Xobs, the contribution of conformally invariant fermions ~ Xobs (TeV)'^, studied in 

Section II, is much more important than that of gravitons, which scales as A^^^. 



We have shown that in brane-world scenarios with a warped extra dimension, it is in 
principle possible to stabilize the radion through the Casimir force induced by bulk fields. 
In particular, conformally invariant fields induce an effective potential of the form ( |3.9| ) as 
measured from the positive tension brane. From the point of view of the negative tension 
brane, this corresponds to an energy density per unit physical volume of the order 



where A is a calculable number (of order 10 ^ per degree of freedom), and A ~ 0/(M'^£)^/^ 
is the dimensionless radion. Here M is the higher-dimensional Planck mass, and £ is the 



V. CONCLUSIONS AND DISCUSSION 



V 



eff '>'[(! -A) 



AX^ 



- + a + pX'^ , 
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AdS radius, which are both assumed to be of the same order, whereas rupi is the lower- 
dimensional Planck mass. In the absence of any fine-tuning, the potential will have an 
extremum at A ~ 1, where the radion may be stabilized (at a mass of order nipi). However, 
this stabilization scenario without fine-tuning would not explain the hierarchy between nipi 
and the TeV. 



A hierarchy can be generated by adjusting (3 according to (|3.11|) , with Xobs ~ 
{TeV/rripi) ~ 10~^^ (of course one must also adjust a in order to have vanishing four- 
dimensional cosmological constant). But with these adjustement, the mass of the radion 
would be very small, of order 

ml ~ Xobs M-H-"" ~ XobsiTeVf. (5.1) 

Therefore, in order to make the model compatible with observations, an alternative mecha- 
nism must be invoked in order to stabilize the radion, giving it a mass of order TeV . 

Goldberger and Wise [§[0, for instance, introduced a field v with suitable classical po- 
tential terms in the bulk and on the branes. In this model, the potential terms on the branes 
are chosen so that the v.e.v. of the field in the positive tension brane v+ is different from the 
v.e.v. on the negative tension brane f_. Thus, there is a competition between the potential 
energy of the scalar field in the bulk and the gradient which is necessary to go from to f _. 
The radion sits at the value where the sum of gradient and potential energies is minimized. 
This mechanism is perhaps somewhat ad hoc, but it has the virtue that a large hierarchy 
and an acceptable radion mass can be achieved without much fine tuning. It is reassuring 
that in this case the Casimir contributions, given by ( ^.1|) , would be very small and would 
not spoil the model. 

We have also calculated the graviton contribution to the radion effective potential. Since 
gravitons are not conformally invariant, the calculation is considerably more involved, and a 
suitable method has been developed for this purpose. The result is that gravitons contribute 
a negative term to the radion mass squared, but this term is even smaller than ( |5.1| ), by an 
extra power of Xobs- 



In an interesting recent paper |]I2[, Fabinger and Horava have considered the Casimir 
force in a brane-world scenario similar to the one discussed in this paper, where the internal 
space is topologically jZ^- In their case, however, the gravitational field of the branes 
is ignored and the extra dimension is not warped. As a result, their effective potential 
is monotonic and stabilization does not occur (at least in the regime where the one loop 
calculation is reliable, just like in the original Kaluza-Klein compactification on a circle 
|3^). The question of gravitational backreaction of the Casimir energy onto the background 



geometry is also discussed in [|T^]. Again, since the gravitational field of the branes is not 
considered, they do not find static solutions. This is in contrast with our case, where static 
solutions can be found by suitable adjustment of the brane tensions. 

Finally, it should be pointed out that the treatment of backreaction (here and in |T^) ap- 
plies to conformally invariant fields but not to gravitons. Gravitons are similar to minimally 
coupled scalar fields, for which it is well known that the Casimir energy density diverges 



near the boundaries 1T3[. Therefore, a physical cut-off related to the brane effective width 
seems to be needed so that the energy density remains finite everywhere. Presumably, our 
conclusions will be unchanged provided that this cut-off length is small compared with the 
interbrane separation, but further investigation of this issue would be interesting. 
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APPENDIX A: EVALUATION OF THE GENERALIZED ZETA FUNCTION 



In this Appendix, we evaluate the generalized zeta function ( |4.26| ) which appears in the 
expression of the effective potential ( |4.25| ). We shall closely follow the method of Ref. pD| , pS| . 
The asymptotic behavior of the function F in the integrand of (|4.26|) is given by 



TlVXt 

—i 



Hence we rewite the zeta function as 
s ' 



(Al) 



cr=± 



27rz I Jc, 



dtt 



In 



-air 



-an 



(A2) 



Here, the contours C± are the upper and lower halves of the contour C in Fig. 3. 

The contribution from the circle at infinity in the first term vanishes, so this part of the 
contour can be dropped. After that, the domain of convergence of the expression is extended 
to 3?s > — 1. The second term can be explicitly evaluated for large s as 



dtt 



-l-s 



In 



/ rftr^-^ln[T7rv^te±(^ 
f /ln[=FivrV^] 



M -s Is, ^(1 1 s 

— \e ' + -e~'\YieT—^ ^£ 

/ s \ — s 



(A3) 



Thus, the whole expression can be analytically continued to the strip given by — 1 < < 0. 
Taking the e — > limit after this analytic continuation, we obtain 

C{s) = ^ sin (^y) dpp-^-' In [nVXpe-^^-^^'Fitp)] 

= ^ sin (^^) dp p-i-^ In [2 v^pe-^i-^)" (p)ir,_i ( Ap) - {Xp)K,_, (p) }] . ( A4) 

The asymptotic expansion of the modified Bessel functions is given by 



^2TTt 



C,{-t)+0{e-') 



(A5) 
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with 

We define coefficients Pr by 



r (z/ + r + i 



r=o r! r ( - r + 2 



fY(2t)-^ 



ln[C,-i(p)] = E, 

r=l P 



Now, we rewrite the expression for ({s) as 



= — sin 

71 



ITS 

Y 



A 



dtt~^-'\n 

POO 

Jo 



2t/7re%_i(t) 



+ / dpp 



In 



27rpe-%_i(p) 

^ _ Vi(Ap) i^,,-i(p) 
ir^_i(Ap) Vi(p) _ 



The first integral in the square brackets can be evaluated in the following way. 
Defining the function 



r=l 



we can rewrite that integral as 



dtt-^-'\n 



2t/Tie^K^_i{t) = J 'dtr^'' In 

roo 

+ dtt 

Jtn 



-1-S 



'to 
D-2 



In 



2t/7i e%_i(t) 



^ s + r J to 



(A6) 



(A7) 



(A8) 



(A9) 



(AlO) 



Except for poles, this expression is now analytic in the strip —{D — 2) < 3?s < 0. The 
integral in the last term is given by 



dtr'-^e-- 



to 



s + {D-l) 



(All) 



After substitution of this expression in (|A10|) , we can perform analytic continuation to 
-{D - 2) > 3?s > -D. Then, taking the to ^ limit, the right hand side of Eq. (^TUp 
becomes 



/•oo 

/ dtr^-' 


(^In 


Jo 





2t/7re%_i(t) 



- m) + 



s + {D-l) ■ 



(A12) 
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Using r(l + x) = 1 — 7a; + 0{x'^), we find 



2t/7re*K^_i(t) 



s + {D-r 



+ Xa' + 0(s + (D-1)), (A13) 



where 



lK = -lf3D-i + l^dtt''-^ [In 



2t/7re%_i(t) 



R{t) 



(A14) 



The second term in the square brackets in ( |A8|) can be evaluated in a similar way as 



dp p-'-^ In [v^e-%_i(p)] = ^ ^^^-^ + X, + o + (/j - 1)) , (A15) 



where 



T/ = -7/9d-i + / dpp 



In 



27rp e^%_i(p) 



i?(-p) . 



Substituting the above results into (|A8|) , we obtain 



1 



ai - D) = -i-iyr^Pn-iil + X 



1-D\ 



+ (D = 3), 

+ (D = 5), 



and 



C'(i-i?) = -(-i)^ 



1 + A 



l-D 



D-l 



A^'^lnA 1 + V(A) 



l-D 



(A16) 



(A17) 



(A18) 



where V(A) is given in ( [4.28|) . Then, substituting these results into ( [4.25| ), we obtain ( [4.27D . 
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